GEOPHYSICAL RESEARCH LETTERS, VOL. ???, XXXX, DOL10.1029/, 



On Zonal Jets in Oceans 

Balasubramanya T. Nadiga 

Los Alamos National Laboratory, MS-B296, Los Alamos, NM-87545 



Balasubramanya T. Nadiga, MS-B296, Los Alamos, NM-87545; balu@lanl.gov 



DRAFT 



February 2, 2008, 12:10pm 



DRAFT 



X - 2 NADIGA: OCEANIC ZONAL JETS 

We find that in parameter regimes relevant to the recently observed alter- 
nating zonal jets in oceans, the formation of these jets can be explained as 
due to an arrest of the turbulent inverse- cascade of energy by free Rossby 
waves (as opposed to Rossby basin modes) and a subsequent redirection of 
that energy into zonal modes. This mechanism, originally studied in the con- 
text of alternating jets in Jovian atmospheres and two dimensional turbu- 
lence in zonally-periodic configurations survives in spite of the presence of 
the meridional boundaries in the oceanic context. 



DRAFT 



February 2, 2008, 12:10pm 



DRAFT 



NADIGA: OCEANIC ZONAL JETS X - 3 

1. Introduction 

A proposed explanation of the alternating zonal jets in Jovian atmospheres is that they 
are due to a tendency of turbulence in thin shells on the surface of a rotating sphere to 
organize itself into zonal jets (e.g., see [Vasavada and Showman, 2005; Galperin et al., 
2004]). The anisotropic jets result from an interplay between an inverse cascade of energy 
[Kraichnan, 1967; Charney, 1947] and the latitudinal variation of the vertical component 
of planetary rotation (e.g., see [Newell, 1969; Rhines, 1975]). While baroclinic instability 
and convective processes are thought to be the main sources of small scale energy, classical 
geostrophic turbulence theory [Charney, 1947] predicts a cascade of this energy (vertically) 
to larger scales as well in a process that has been termed barotropization. Hence, in 
the context of this explanation of atmospheric-zonal jets, they have been simulated and 
studied extensively using the barotropic vorticity equation on either the doubly-periodic or 
zonally-periodic beta-plane or on the surface of a sphere using forced-dissipative settings. 
In these settings, the effect of geometry on dynamics is minimized in the sense that the 
flow in the zonal direction, the direction in which the dynamics of the Rossby waves are 
highly asymmetric, is homogeneous. Dynamically, the formation of the zonal jets in this 
homogenous setting is thought to involve certain kinds of resonant interactions (sideband 
triad and quartet) of Rossby waves packets whose amplitudes are slowly varying functions 
of space and time [Newell, 1969]. 

More recently, observational [Maximenko et al., 2005] and computational [Nakano et 
al., 2005] evidence point to the occurence of multiple alternating zonal jets in the world 
oceans as well. However, the dynamics underlying their formation is not clear. 
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On the one hand, given that the governing equations are the same in the atmospheric 
and oceanic contexts, it would not be unreasonable to expect, from a turbulence point of 
view, that the same dynamical mechanism — Rossby wave dispersion arresting the inverse 
cascade of energy and redirecting it into zonal modes — underlies the phenomenon, be 
it in the ocean or in the atmosphere. Clearly, unlike the constant stratification of the 
atmosphere, surface-intensified stratification in the oceans inhibits full barotropization 
(e.g., see [Fu and Flierl, 1980]). Nevertheless, the importance of the barotropic mode 
(with a thermocline depth of 1 km in a 5 km deep ocean) is clearly borne out in Fig. 2 
and Table 1 of Fu & Flierl [1980] and other such studies confirm an inverse cascade of 
barotropic kinetic energy. High vertical coherence of jet structure in models [Nakano et 
al., 2005; Maximenko et al., 2005] further suggests the importance of barotropic dynamics. 

On the other hand, the presence of boundaries can, besides being able to support viscous 
boundary layers and act as sources/sinks of enstrophy, allow for new (inviscid) mecha- 
nisms. For example, in a closed basin, (a) Fofonoff gyres arise as statistical equilibrium 
solutions of the barotropic vorticity equation, and (b) Rossby basin modes arise, reso- 
nant interactions of which have been studied as mechanisms for generating both mean 
flows (starting with Pedlosky [1965]) and mesoscale variability (starting with [Harrison 
and Robinson [1979]). Such mechanisms could possiblly generate alternating zonal jets 
as well. Interestingly, LaCasce [2002] finds that the arrest of the inverse cascade of en- 
ergy by basin normal modes is largely isotropic. However, in a recent article studying 
rectification processes in a three layer quasi-geostropic beta plane basin, Berloff [2005] 
concludes that the alternating zonal jets he found in that setting were most likely driven 
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by nonlinear interactions between some meridionally structured baroclinic basin modes 
and some secondary (i.e., related to finite amplitude background flows) basin modes. If 
this were to be the most important mechanism for the formation of alternating zonal jets 
in ocean basins, by involving spatially-nonlocal (basin) modes this mechanism would be 
fundamentally different from the (spatially) local arguments of turbulence that are usually 
thought to apply in the atmospheric context. 

In this letter, we demonstrate that in parameter regimes relevant to alternating zonal 
jets in the oceans, such jets can be formed by free Rossby waves (as opposed to Rossby 
basin modes) arresting the inverse-cascade of energy. We then go on to show that the 
jet width scales well with Rhines' scale. This suggests that the dynamics of alternating 
zonal jets in oceans are likely local and in this sense similar to those in previously studied 
atmospheric contexts. We suggest that the nonlocal resonant-interaction-of-basin-modes 
mechanism becomes more important at larger values of turbulent kinetic energy (TKE). 
Curiously, only the latter regime has been investigated before within the framework of 
the barotropic vorticity equation [LaCasce, 2002], and as far as we know this is the first 
time that alternating zonal jets have been obtained in a closed basin using the barotropic 
vorticity equation. 

The rest of the letter is structured as follows: The next section briefly describes the 
modeling approach, and the following one presents computational results. As a matter 
of convenience, and with no loss of generality, these two sections consider the governing 
equations and present the results in a nondimensional form. The final section establishes 
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the correspondence between the nondimensional parameter values considered and their 
dimensional counterparts in actual ocean settings. 



2. Modeling Approach 

We consider the barotropic vorticity equation 

f t +J(^ q ) = F + D (2.1) 

for the evolution of barotropic potential vorticity q — £ + fly — V 2 ip + fly, where ( is 
relative vorticity, ip is velocity streamfunction, F is forcing, D is dissipation and J(, ) is the 
Jacobian operator given by J(ip, q) = — f^ff + §^f§- The above equation is considered on 
a midlatitude beta plane with a latitudinal gradient of the vertical component of rotation 
of fl\ y-coordinate increases northwards and x-coordinate eastwards in a closed square 
basin, 2n on a side, discretized into 1024x1024 cells. An energy and enstrophy conserving 
finite-differencing is used with Runge-Kutta timestepping [Greatbatch and Nadiga, 2000] . 
Given the inverse-cascade of energy of 2D turbulence, forcing F is concentrated around a 
high wavenumber kf, as a combination of sines and cosines consistent with the boundary 
conditions used. Their amplitudes a, drawn randomly from a Gaussian distribution, are 
delta-correlated in time resulting in F = a(t)/ V5tf(kf,t) with an energy input rate e of 
° 2 IS fV ^fdxdy and an enstropy input rate r] of kje. In all the computations presented, 
given the domain size of 27rx27r discretized into 1024x1024 cells, k max is 512 and kf is 
between 128 and 129. 

Dissipation D = — z/ p V 2p C — ^oC> consisting of a small-scale-selective component to 
dissipate the (largely) downscale-cascading enstrophy input at the forcing scale, and 
Rayleigh friction component that mainly acts to dissipate the (largely) upscale-cascading 
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energy. At lateral boundaries, besides no through-flow, we use superslip boundary condi- 
tions. The coefficient u p is diagnosed dynamically in terms of the enstrophy input rate as 



Given the above setup, the problem consists of three important parameters: (3, e and 
vq. We briefly recall a few relevant spatial scales in terms of these parameters. First, in 
purely two-dimensional turbulence, a Kolmogorov scale for the dissipation of energy may 
be obtained using the usual arguments as 



(e.g., see [Danilov and Gurarie, 2002]). In the absence of (3 this would be the scale at which 
Rayleigh friction would act to stop the inverse cascade of energy. However, in the presence 
of P, Rossby wave dispersion can instead arrest the inverse cascade and redirect energy 
into zonal modes. In the absence of large scale friction, and under the assumption that the 
spectral flux of energy in the inverse-cascade inertial range is determined by the energy 
input rate e (due to forcing), this would happen at kp = (/^/e) 1 / 5 [Vallis and Maltrud, 
1993]. If, however, energy is concentrated near kp, this arrest mechanism would occur at 
the Rhines' scale kp = \J f3/U rms [Rhines, 1975] (also obtained by equating the turbulence 
frequency U\kp\ and the Rossby wave frequency — /3cos<f>/\kp\, where = tan -1 k y /k x ). 
Given the largely upscale-cascading nature of energy, the small-scale-selective dissipation 
operator plays a relatively minor role in dissipating energy, so that dE/dt ~ e — 2u E, 
with energy leveling off at about e/2z/ . Using this energy balance in the expression for 
Rhines' scale leads to [Danilov and Gurarie, 2002; Smith et al., 2002] 



v p = CkV 3 Arc 2p ', using Kolmogorov-like ideas and a Kolmogorov scale of Ax. 




(2) 




(3) 
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3. Results 

Table 1 gives the basic parameters and the derived scales discussed above for a series 
of simulations. In all the cases considered, care is taken to ensure that the spectrum of 
the zonal component of energy has equilibrated. While there are important differences 
between some of the cases in Table 1, we postpone a detailed discussion of these differ- 
ences to a later article, and go on to examine a representative case — case C presently. An 
examination of the instantaneous, zonally-averaged, zonal-velocity and relative-vorticity 
fields plotted as a function of latitude in Fig. 1 suggests alternating zonal jets of a charac- 
teristic width. To further verify this, we examine a few other familiar diagnostics. First, 
Fig. 2 shows the time-mean two-dimensional zonal-velocity field after the flow has reached 
statistically-stationarity, and the alternating zonal jets are evident in this figure. Note that 
a) even though the forcing is homogeneous, the jets are more pronounced to the west, b) 
unlike with observations, time-mean jet signatures are obtained and analysed, and c) the 
geometry of the jets are not significantly different when the time- varying components are 
analysed (not shown). Finally, while the meridional gradient of time-averaged potential- 
vorticity is dominated by (3 (stable), the instantaneous flow quite frequently violates the 
barotropic stability criterion u yy < (3. That these alternating zonal jets are related to 
anisotropization of the inverse cascade of energy of two dimensional turbulence by Rossby 
wave dispersion is verified by the the dumbell shape near the origin, characteristic of the 
process (e.g., see [Vallis and Maltrud, 1993]), in the contour plot of the two dimensional 
spectral density of energy in Fig. 3. 
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It is not our intent to verify various universal scalings of spectra in this problem, but to 

use it as a diagnostic to further confirm the nature of the dynamics. To this end, we show 
in Fig. 4 the range of spectra that we obtain in the parameter range considered. These 
figures show the ID energy spectra averaged over an angle of n/6 around = (residual 
flow) and = n/2 (zonal flows) [Chekhlov et al., 1996]. Both the residual and zonal 
spectra have futher been compensated for the A;~ 5 / 3 scaling. (A compensation for e 2 / 3 — 
following Kolmogorov scaling E(k) = Ck€ 2 / 3 k~ 5 / 3 — is avoided since while that would be 
appropriate for the residual component, it would not be appropriate for a possibly different 
scaling of the zonal component such as E z (k) = C z f3 2 k~ 5 . However, e 2 / 3 compensation has 
been applied to the residual spectra to establish the value of the Kolmogorov constant 
Cfc). A common and important feature of all the cases is that in the inverse-cascade 
regime, while at the high-wavenumber end, the zonal and residual spectra scale similarly, 
at lower wavenumbers the zonal spectra lie above the residual spectra and show steepening 
before they peak. This behavior is as expected and verified by various investigators in 
the periodic case relevant to the atmosphere. Furthermore, like in the computations of 
Danilov and Gurarie [2001], our spectra display significant non-universal behavior. For 
example, while in cases A and F, the residual spectra clearly verify the classic Kolmogorov 
scaling C k e 2 / 3 k~ 5 / 3 with a Kolmogorov constant C k of about 6, that is not the case for 
cases C and J. We note parenthetically that (a) the distribution of spectral energy flux 
as a function of wavenumber (not shown) bears remarkable resemblance to that derived 
using Aviso, TOPEX/Poseidon, and ERS-1/2 data [Scott and Wang, 2005], and (b) on 
using the spectral flux of energy as a function of wavenumber, as opposed to a constant 
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value, Ck remains close to 6 in the high wavenumber range of the inverse- cascade, but 
then begins to rise at the lower wavenumbers. 

Next, we identify the jet width with the wavenumber at which the (uncompensated) 
zonal-spectrum peaks, k p , and verify it by referring to physical-space pictures like in 
Figs. 1 and 2. This number is recorded for each of the cases in the last column of Table 1. 
We note, that the wavenumber at which the residual spectrum peaks is close to this 
wavenumber as well. In Fig. 5, we plot the above measure of jet width (k p ) against the 
Rhines' scale (kp) and find excellent agreement, like in the periodic (atmospheric) case 
(e.g., see Vallis and Maltrud [1993] and Danilov and Gurarie [2001]. 

4. Discussion 

The simulations and analyses presented in the previous section show clearly that there 
are parameter regimes wherein the dynamics of the alternating zonal jets in a midlatitude 
ocean basin are not controlled in a fundamental manner by meridional boundaries. That 
is to say, in these parameter regimes, the dynamics of the jets are governed largely by 
spatially local interactions, and the arrest of the inverse-cascade is mediated by free Rossby 
waves as opposed to Rossby basin modes. However, we still need to check if such a 
parameter regime is of relevance to the oceans in order to establish the importance of this 
mechanism in the oceans. 

The pronounced signature of the observed jets in western regions of ocean basins (Fig. 
1 of [Maximenko et al., 2005]) would generally be attributed to the elevated levels of TKE 
in the separated western boundary current (WBC) regions. However, our simulations use 
homogeneous forcing but still display similar enhanced jet signatures in the west. This 
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leads us to suspect that the enhanced signature of the jets in the west is more due to its 
internal dynamics, rather than due to the elevated levels of TKE in the WBC regions, 
and that the jets are controlled more by the ambient (lower) levels of interior TKE. 
An approximate range of 25 to 100 cm 2 /s 2 is obtained for the latter by examining an 
altimetry-derived TKE map for the North Atlantic [pers. comm.; Scott, 2006]. Keeping 
this in mind, first consider case C discussed extensively above: the peak wavenumber k p is 
26 (Table 1); using the observed [Maximenko et al., 2005] dominant wavelength of 280 km. 
leads to L re f of 1160 km. Then, using a typical midlatitude value of (3 re f of 2 10 -11 m _1 s _1 
and an r.m.s. value of 7.5 cm/s (mid-range) for the domain-averaged interior geostrophic 
velocity anomaly (TKE), leads to a f3 n d [=(5 re fL 2 re ^ /U re f) of 360, corresponding well with 
320 used for case C. As for the ranges of parameters considererd, 5 < k p < 55 (Table 
1); using a range of wavelengths of 500 to 250 km and do main- averaged interior TKE 
level corresponding to 5 < U rms < 10 cm/s, gives (3 nd in the range 30-1900 (c.f. range 
of 80-1280 in Table 1). In light of this, we suggest that the local mechanism wherein 
the arrest of the inverse cascade is mediated by free Rossby waves may be important in 
explaining the formation of alternating zonal jets in the world oceans. Obviously, further 
work is necessary to definitively establish the relevance of this mechanism to the oceans. 
The author thanks Boris Galperin for discussions and the referees for their criticism. 
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Table 1. Basic parameters, derived scales and the jet-width wavenumber for the ten simulations 



considered 



Case 




e 






kfr 




kp 



A 


80 


0.50 


0.1 


15.9 


2.2 


15.9 


17.5 


B 


160 


0.50 


0.1 


24.1 


2.2 


22.2 


15.5 


C 


320 


0.50 


0.1 


36.6 


2.2 


31.3 


26 


D 


640 


0.50 


0.1 


55.5 


2.2 


38.5 


37 


E 


1280 


0.50 


0.1 


84.0 


2.2 


66.7 


55 


F 


80 


33.5 


0.1 


6.0 


0.2 


5.96 


4.5 


G 


80 


128 


0.4 


4.6 


0.8 


5.88 


7 


H 


1280 


126 


0.4 


27.8 


1.1 


22.2 


20 


I 


1280 


296 


0.4 


23.4 


0.7 


18.2 


16 


J 


1280 


513 


0.4 


21.0 


0.6 


15.9 


13 
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Figure 1. Meridional plot of the instantaneous zonally- averaged zonal- velocity and vorticity. 
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Figure 2. The alternating zonal jets are evident in the time-averaged, two-dimensional zonal- 
velocity field. While forcing is homogeneous, jets are more prominent in western regions. 




Figure 3. The time-averaged two-dimensional energy spectrum displays the familiar anisotropic 
'dumbell' shape. 
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Figure 4. Time-averaged one-dimensional zonal (dot-dashed line) and residual (solid line) 
energy spectra. Both have a k~ 5 ^ 3 compensation. See text for details. 
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Figure 5. A plot of the wavenumber at which the zonal-spectrum peaks, k p against the Rhines' 
scale k§. Symbols correspond to the cases in Table 1 and line to the linear least squares fit. 
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